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We consider cosmological evolution in Brans-Dicke theory with a population of primordial black 
holes. Hawking radiation from the primordial black holes impacts various astrophysical processes 
during the evolution of the Universe. The accretion of radiation by the black holes in the radiation 
dominated era may be effective in imparting them a longer lifetime. We present a detailed study of 
how this affects various standard astrophysical constraints coming from the evaporation of primor- 
dial black holes. We analyze constraints arising from the present density of the Universe, the present 
photon spectrum, the distortion of the cosmic microwave background spectrum and also from pro- 
cesses affecting light element abundances after nucleosynthesis. We find that the constraints on the 
initial primordial black hole mass fractions are tightened with increased accretion efficiency. 
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I. INTRODUCTION 

The Standard Model of Cosmology invokes General Theory of Relativity(GTR) as the theory of gravity. 
However, all the confirmative tests of GTR have been carried out at low energy. This has led people to 
believe and explore the deviations from GTR at high energy regimes like cosmic evolution at very early 
times. Brans-Dicke (BD) theory of gravitation [l[ stands out as one of the most attractive alternatives 
to GTR because it involves minimal extension over GTR through introduction of a scalar field <f>. In the 
BD theory the gravitational constant becomes function of time and is proportional to the inverse of the 
scalar field <f> which couples to gravity with a coupling parameter lo. The more general nature of BD 
theory is evident from the fact that in the limit u — > oo it goes over to GTR. Solar system observations 
require uj > 10 4 Q. The ubiquitous nature of BD theory is also evident from the fact that it appears in 
the low energy limit of Kaluza- Klein and String theories Thus, BD theory has been used for tackling 
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a number of cosmological problems such as inflation [4( , early and late time behaviour of the Universe |5[ , 
cosmic acceleration and structure formation Q, coincidence problem Q and problems relating to black 
holes 0. 

It has been pointed out that in the early stage of the Universe Primordial Black Holes (PBHs) could be 
formed due to various mechanisms such as inflation Q , initial inhomogeneous conditions [ltj EH , phase 
transition 12], bubble collision [l3| and decay of cosmic loops [l4j]. In the usual formation scenarios, 
the typical mass of PBHs at the formation time could be as large as the mass contained in the Hubble 
volume Mh ranging down to about 1CU 4 Mh [Hj]. The formation masses of PBHs could be thus small 
enough to have evaporated completely by the present epoch due to Hawking evaporation [l6|. However, 
since the cosmological environment is very hot and dense in the radiation-dominated era, it is expected 
that appreciable absorption of the energy-matter from the surroundings could take place. It has been 
noticed that such accretion is most effective in altered gravity scenarios where the PBHs grow due to 
accretion of radiation at a rate smaller than that of the Hubble volume, thus providing for enough energy 
density for the PBHs to accrete causally. This is responsible for the prolongation of the lifetime of PBHs 
in braneworld models !17] as well as in scalar-tensor models [l8j . 

The feasibility of black hole solutions in BD theory was first discussed by Hawking [l9| . Using scalar- 
tensor gravity theories Barrow and Carr [2(| have studied PBH evaporation during various eras. It 
has been recently observed [2l| that in the context of Brans-Dicke theory, inclusion of accretion leads 
to the prolongation of PBH lifetime. Once formed, these PBHs will influence later cosmological epochs, 
leading to a number of observational constraints on their allowed abundance. These have been extensively 
investigated in the case of standard cosmology [22|-[24[ . In a recent work, Carr et al. [25[ have performed 
a detailed numerical study of the effect of the emission of quark and gluons by PBHs on the standard 
constraints. The standard constraints also get altered in different gravitational theories, as was studied 
in the context of Brans-Dicke cosmology without accretion [2(J, and brane- world cosmology [261-2 8]. The 
aim of the present paper is to reanalyse the main constraints in the Brans-Dicke theory with the inclusion 
of accretion of radiation in the early Universe. The standard constraint formalism could get modified due 
to the change in theory of gravity and also due to the effect of accretion at early times, as was shown in 
the case of braneworld gravity by Clancy et al. [26j j . In the present study we extend our previous work on 
the evolution of PBHs in BD theory including accretion [21j to obtain several astrophysical constraints 
following the style of Clancy et al. |26j . 

The plan of the paper is as follows. In the next section we first provide the key expressions related 
to PBH evaporation, accretion and lifetimes in the BD scenario, and then discuss the modified con- 
straint formalism in BD theory that we use subsequently. In Section III we consider the observational 
constraints obtained from the present density of the Universe. In Section IV we discuss the present 
photon spectrum and constraints following from it. The constraints arising out of the distortion of the 
Cosmic Microwave Background (CMB) spectrum are evaluated in Section V. Light element abundances 
and photo-disintegration of the deuterium nuclei lead to further constraints on the initial PBH mass frac- 
tion, which are presented in Section VI. A summary of our results along with a Table with quantitative 
estimates of the various constraints are presented in Section VII. 



II. THE CONSTRAINT FORMALISM IN BRANS-DICKE THEORY 



For discussing the constraints that can arise from PBH evaporation, we label an epoch by cosmic time 
t. PBHs which are not evaporated by time t will only contribute to the overall energy density of PBHs. 
As the present observable Universe is nearly flat and, therefore, possesses critical density, the PBH mass 
density can be constrained on the ground that it should not overdominate the Universe. PBHs evaporate 
by producing bursts of evaporation products. The Hawking radiation from the PBHs which evaporate well 
before photon decoupling will thermalize with the surroundings, boosting the photon-to-baryon ratio (29j . 
In the case of evaporation after photon decoupling, the radiation spectrum is affected and subsequently 
redshifts in a monotonic manner. Thus, constraints arise from the cosmic background radiation at high 
frequencies [30|, HH EJ] ■ Further, if the PBHs evaporate close to the time of photon decoupling, it cannot 
be fully thermalised and will produce distortion in the cosmic microwave background spectrum. Generally 
speaking, at a given epoch, the constraint on various physical observables is usually dominated by those 
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PBHs with a lifetime of order of the epoch in question. Hence, the observational constraint can be 
translated into an upper limit on the initial mass fraction of PBHs. 

We consider a flat FRW Universe which is radiation dominated upto a time t e and matter-dominated 
thereafter. In the radiation-dominated era [2(3] 

«(i)«i 1/2 ; G = G (^)", (1) 

where Go denotes the present value of G, and n is related to the Brans-Dicke parameter w by n = 
2/ (4 + 3w). [In view of the observational bound on w 0, we consider the value of n as n ~ 0.0001 in our 
subsequent calculations] . The radius and temperature of a PBH is given by 

r = 2G m n M; T BH = . (2) 

Similarly, in the matter-dominated era 

a(t)oc^ 2 -")/ 3 ; G(i)=G (y) n , (3) 

and 

In Brans-Dicke theory the evolution of a PBH [21| is described by 

dM _ fdM\ rdM\ 
dt V dt J (ICC \ dt J evap 

with the accretion and evaporation rates given by 



and 



where a = er/(2567r 3 Gg), with cr the Stefan-Boltzmann constant. / is the accretion efficiency which 
depends primarily on two factors: (i) how small the size of the black hole is compared to the cosmological 
horizon, and (ii) how efficient the black hole is in absorbing background radiation. The first factor 
is contingent upon various formation mechanisms. The second factor could depend, in turn on the 
thermal properties of the surrounding radiation, the mean free path of absorbed particles, as well as the 
backreaction due the black hole on the cosmological metric. Due to several uncertainties in the estimation 
of the above quantities, it is usually accepted to consider an accretion efficiency / in the range between 
and 1. 

Accretion makes the mass of PBH to grow a maximum value [2lj 

M max = M{t c ) = — -4-t (8) 

where Mi is the initial mass of PBH and t c is the time at which evaporation rate is equal to the accretion 
rate. Thus / is further restricted to lie in between and | in Brans-Dicke theory. We assume the 
standard mechanism of PBH formation due to the gravitational collapse of density perturbations at the 
cosmological horizon scale [jjj . which leads to Mi ~ G~ 1 (ti)ti. 
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The evaporation time of PBHs which completely evaporate in the radiation dominated era, is given by 

m 



The PBHs which evaporate in the matter dominated era have their initial mass and evaporation time 
related as [H| 

r 3 •> r // \2n r r // \2n+l n 1/3 

"i-f'-i'l'Hs) ..{1+P.+D- 1 • < 10 > 

Constraints on the allowed abundance of PBHs of a certain lifetime are formulated as upper bounds 
on their mass fraction. This mass fraction at (Mi) is defined as the ratio of the energy density due to 
PBHs of initial mass Mi and the background radiation density, at a time t >U, as 

a t (M i ) = j-—. (11) 

Prad(t) 

The black holes first gain mass through accretion up to time t c . During this stage at scales as (Ma). 
Subsequently, they begin to loose mass very gradually (i.e., the mass stays nearly constant, and therefore 
at scales as a). This stage comes to an end when at reaches the value a evap beyond which the black 
hole looses most of its mass in a final burst of evaporation. It then follows that if initial and final mass 
fractions denoted by and a evap respectively, are related by 



M(t c ) a(t eva p ) 
l M(U) a(U 



%evap — a% n j \ n f± \ ' V / 



Using equation (j8]), we can write 
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a(tevap) 



The purpose of the following sections is to reconsider observational constraints on a evap at different 
cosmological epochs, trace them back to obtain constraints on the initial mass fraction 014 taking into 
consideration accretion of radiation in the early universe by the PBHs in Brans-Dicke theory. It may 
be mentioned here that the recent paper by Carr et al. (25| has reanalyzed the standard constraints 
on primordial black holes by considering the effects of emission of quarks and gluons and the resultant 
secondary emission of photons. It was shown that the effect of secondary photon emission could alter the 
standard constraints on PBH fraction by a couple of orders of magnitude in certain cases, e.g., deuterium 
constraint, while leaving the standard constraints more or less unaltered in certain other cases, e.g., 
distortion of CMB spectrum. The results obtained by Carr et al. j25[ are based on detailed numerical 
analysis. Of course, such a scenario of emission would also impact constraints on Brans-Dicke primordial 
black holes in more or less similar ways as they impact PBHs in the standard picture. However, the lack of 
analytical results describing the effect of such emission on the constraint formalism makes it considerably 
harder to perform a similar analysis in the context of an altered gravitational scenario. Since the primary 
aim of the present paper is to study the effect of Brans-Dicke theory on the constraint formalism, we do 
not attempt to perform any similar numerical analysis here. In stead, we make a note of the effects that 
emission of secondary photons could have on the different constraints, based on the results of Ref. [25j . 



III. THE PRESENT MATTER DENSITY OF THE UNIVERSE 



For a particular value of / in Eq.(j8]) with t evap — to, and using Eq.lfTOf one gets the formation time of 
PBHs that are evaporating today. Black holes formed later are essentially still intact and their density 
is constrained by the observed matter density in the present Universe. Here we are tracking the relative 
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densities of PBHs to radiation and we must ensure that, given the observed radiation density, this ratio 
does not imply that the PBH density exceeds the observed matter density of about 0.3 of the critical 
density. Phrased in this way, the constraint applies regardless of the presence of a cosmological constant 
and indicates that for any PBHs surviving to the present we must have 

«o(M) < ™ (14) 

"7,0 

The cosmic microwave background corresponds to a photon density of f2 7i oft. 2 = 2.47 x 10 -5 (with h = 0.7) 
and conservatively, we can ignore the cosmic neutrinos. Thus, for PBHs that are about to evaporate today, 
tevap > <o, we get 

3 

a eV ap < pr— ~ 6 x 10 3 . (15) 

"7,0 



a, <3.43x 10- 18 x (l-|/) 3/2 . (16) 



Using eauation (|13[) . one can find 

rv, < 3 43 x 10~ 18 x (\ - 

The constraint on the initial PBH mass fraction in the tandard cosmology is obtained as [3l[ 

oti < 10" 18 . (17) 

IV. THE PRESENT PHOTON SPECTRUM 

If PBHs evaporate between the time of photon decoupling^ec) and the present day, their radiation 
spectra will not be appreciably influenced by the background Universe, apart from being redshifted. Thus 
the spectra could constitute a fraction of the cosmic background radiation. The number of particles of 
a certain species j, emitted in 4-dimensional spacetime by a black hole of temperature Tbh, in a time 
interval dt and with momentum in the interval (k, k+dk) is 

, _ « dt d k , . 

^'"'"^M (18) 

with uj 2 — k 2 + m 2 , where Cj(fc) is the emission cross-section for species j with momentum k and m is 
the mass of the particle. 

For emission of photons with energy in the interval (E, E+dE) above equation becomes 

^-^ ^X)-' " 8 * (19) 

Thus the spectral photon number emitted by a small black hole with lifetime t evap is obtained from 

dN f u " a » a(E) E 2 

But, in the high frequency limit E » Tbh, all cross section reduce to the same value [26| as 

1 

Now, the spectral number becomes 



- = ^ = 4^- (21) 



_ = _ / KTbh ' dt (22) 

dE 8^J esp(_|L)_i ■ 1 > 



G 



In the high energy limit, the above integral yields the spectrum [26| 

dN 24.886 , ft e \ n fi „ , 

(where a has been defined below Eq.Q.) The spectrum declines as E~ 3 like the standard spectrum. 

Now, consider black holes evaporating at a time t evap > td ec - We make the approximation that all the 
energy gets released instantly, but take the spectrum into account. The black hole mass fraction just 
before evaporation is given by 



M(t c ) a(t evap ) 
'—-^ Ml a{U) 



while the number density (which scales as a~ 3 ) is 



, n «(*c) 3 



n e vap(Mi) = n eV ap(M c ) _ u °^ , (25) 



where 



leading to 



fH/T\ Prad(tevap) ma\ 

n evap (M c ) = a evap — — - — (26) 



/ \ Prad{tevap) (- 3 \ 3 / 2 . . 

npBH(tevap) = Oievap Jj7jT\ ~ 2 ) ' 

The energy density in photons of energy E emitted between t evap and t evap + dt evap is [26j j 



dUtev(E) — npBH(t evap )E 2 ——(t evap ) evap , (28) 



We require the present total energy density in Hawking photons at a certain energy scale -Eoi denoted as 

'"(,(/•.,;• 

The photon energy E emitted at time t undergoes redshift due to expansion of the Universe. So the 
photons emitted at time t evap with energy E is related with presently observed energy as 

E(E ) = E J^L (29) 



Since 



one gets 



Prad(tevap) = Prad(t())— 77 (30) 



U (E )= I dUo(E a )= I \™ a *l dU tev [E(E )}. (31) 



Substituting the relevant expressions and using Eq. (fT0|) . we get 



[/»(£„) ="v(£) ~ (( t )* + "(^)'(3o)-"V™J«»)£„ 2 fr" )/3 (2"+l)' 
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The number spectrum of a black hole of initial temperature Tbh peaks at an energy E — bTsH 



with b s=s 5 |33| in the standard treatment. In the BD case the b- value remains unchanged since in 
its determination from the power spectrum of the emitted particles, G which in this case is ~ gets 
cancelled. Therefore, unless a, is sharply peaked at particular initial epochs, the main contribution to the 
integral in Eq.p2p is obtained when E(Eq) = 6Tb# (i e „ ap ), i.e., from PBHs evaporating at t evap = t ma i n , 
where 



tr. 



f En \3/l-3n 



The contribution from PBHs evaporating earlier will come from the high frequency end of their spectrum, 
while PBHs evaporating at later times will contribute radiation that originated in the low frequency end. 
Using Xi = b « 5, the number spectrum equation (|22[) becomes 

dN 14.3 m e pl ft e y 3 /a(t evap )\3 



- 409^6 — (H E o • ( 34 ) 



dE(E ) yevapj 4096tt 6 a \t a J V a(t ) 
Using this equation in equation (|32j) . we estimate the total energy density at energy Eo as 

o( o) " 4096^6 a \T) mpl \TJ (*e) a+n (3a) 1/2 p ra d(to)E H > , (2n+l)-2 

" "dec ~ ' 

which gives 

£/ (#o) oc4 3 - 16n)/(2 - 6n) . (36) 

In the standard spectrum Uq(Eq) varies as E^ 2 . 

Radiation at lower frequencies will originate from the low frequency ends of the instantaneous spectra, 
with the dominant contribution coming from PBHs evaporating around td ec - Its intensity can generically 
be neglected as compared to the main frequency range [26j. For energies En > &Tb#(£o), the dominant 
part comes from the high frequency tail of PBHs evaporating today. The number spectrum equation in 
the high energy limit is used to obtain 

U (E ) = 6.16 x 10 54 ^(7V 4 (l - |/) 3 keVcm- 3 s- 2 keV- 1 . (37) 

The spectral surface brightness I (En), an observational quantity, is related to the integrated energy 
density U (E ) by ^ 

The overall peak in the present spectrum is at E pea k — 6Tbh(*o) ■ So equation (l3"8l gives 

I{E peak ) = 2.08 x 10 20 (l - ^/) 1 a i keVcm~ 2 s' 1 sr^keV' 1 . (39) 

Considering a value of E pea k ~ 100 Mev [3J], which is not too far from the present observational range 
[35l ] of 7-rays, we have 

I obs = 1.11 X 10" 5 keVcm^s^sr^keV- 1 . (40) 
The constraint I(E peak ) < I b s then results in an upper limit on the initial mass fraction as 

1-26 ( - ^ 



a,<5.34xl0-*x 1--/ . (41) 
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For comparison, the corresponding constraint in the standard case is obtained from the gamma-ray 
background at -E peQ & = lOOMeV and reads a, < 10 -27 [23|, [34|, |4l[ . It may be noted here that the recent 
work of Carr et al. [25[ confirms earlier results [36[ that the spectrum of secondary photons is peaked at 
E w 68MeV, independent of the PBH temperature. Our constraint on the PBH mass fraction originates 
from considering the peak value E s» lOOMeV for primary photons. As shown by Carr et al. (25j . 
secondary photon emission may dominate when the PBH mass falls below the QCD scale, in which case 
the constraints may be altered at most by an order of magnitude. Similar considerations would apply to 
the Brans-Dicke case, but the exact magnitude of the constraint would require numerical analysis to be 
evaluated. 



V. DISTORTION OF THE COSMIC MICROWAVE BACKGROUND SPECTRUM 

Hawking radiation emitted at redshifts z < 2 x 10 6 or t > 4 X 10~ 10 to can not be fully thermalised and 
will disturb the Cosmic Microwave Background spectrum. The corresponding modification of the Planck 
spectrum is described by the chemical potential /x (37j . which is related with injected energy density as 

0.71m, (42) 

Prad(t) 

where p evap is the energy density injected by evaporating PBHs. Observational results [38j suggest an 
upper limit on /x given by 

H < 9 x 10~ 5 . (43) 

Assuming that about a half of the energy is emitted in the form of particles capable to disturb the CMB 
spectrum, one can write 

~<W P = 0.71/i. (44) 

Using equation (|4*3"|) . we get 

a evap < 1.28 x 10" 4 , (45) 
which leads to the initial PBH mass fraction constraint as 

3 A 3/2 



m, .. 1.28 ■ 10 -21 x ( 1 } 



/) ■ (46) 



In the standard cosmology the corresponding constraint is [391 ] 

oti < 10~ 21 . (47) 

In a recent work it was shown by Tashiro and Sugiyama (40j , that secondary photon emissions leading 
to non-zero chemical potential for photons could impact the CMB spectrum. However, it follows from 
the analysis of Carr et al. [25|, that the constraint on the PBH mass fraction stays around the value 

io- 21 . 



VI. NUCLEOSYNTHESIS CONSTRAINTS 



Standard big-bang primordial nucleosynthesis is one of the most well-understood processes in the early 
Universe. Therefore, this era is an important benchmark to look for effects due to the interactions of 
particles emitted by PBHs 41]. Several detailed investigations of PBHs in standard cosmology have 
computed the predicted changes in the density of light elements [13, E3] • Existing observational limits on 
the light element abundances have then been used to put constraints on the size of such modifications, 
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which in turn lead to constraints on the numbers of PBHs that could evaporate both during and after 
nucleosynthesis. In the standard scenario, the PBHs which evaporate during nucleosynthesis with t evap 
between Is to 400s, have initial masses in between 10 9 g and 10 10 g. This remains nearly the same in Brans- 
Dicke theory where the initial mass varies between 3.1 x 10 9 x (1 - §/)g and 2.28 x 10 10 x (1 - |/)g. 
Following again the analysis of Clancy, Guedens and Liddle [26| (performed in the context of braneworld 
black holes) here we examine two nucleosynthesis constraints in the context of the Brans-Dickc theory, 
namely the constraint on the increase in production of helium-4 due to the injection of PBH hardons 
[43|,[44[ and the constraint on the destruction of primordial deuterium by PBH photons [43|,|45j]. Note here 
that the nucleosynthesis constraints are the ones that could be most affected by taking into account quark 
and gluon emission by PBHs. Inter-conversion between protons and neutrons due to emitted mesons and 
anti-nucleons increases the n/p freeze-out ratio as well as the final He-4 abundance, as shown by Carr et 
al. [25j in their recent work. Our main aim here is to get a reasonable estimate of how such constraints 
are modified in BD theory. 



A. The Helium abundance constraint 



The total number density of emitted particles from the complete evaporation of PBHs of some initial 
mass may be expressed as 

N em = (48) 

where < E ern > is the average energy of the emitted particles. The ratio of the energy density in PBHs 
at evaporation to the background radiation energy density is therefore 

P pBH - a - < E em > N em 

— <^evap — j-, AT , y^yj 

Prad < i^rad > ^rad 

where < E ra d > and N ra d are the average energy and number density of the particles comprising the 
background cosmological radiation fluid. The ratio of the average energies can be approximated by the 
ratio of the PBH temperature at the onset of evaporation to the background temperature at evaporation, 
i.e., 

T~F ^ = T ■ ( 50) 

\ ^rad ^ evap 

Using equation (2)-(8) and applying the standard cosmological temperature-time relation [46| 

t = 0.3015-^ (51) 

(where g* is a constant having value 10.78 ), one gets 

< E rad > _ T evap , 00w ^ 3 A (M(t c )\-V* 



= 7.28 x 



< E em > Tbh ^ 2 / V m p i 

The total emitted number density at time t evap during nucleosynthesis thus becomes 

:i , ( Mijr) . 2 
mpi 

Using observational estimations, Clancy et. a?.[26| found that 



AT em = 7.28(l-|/)a ewp (^^) Nrad- (53) 



2 8 

Nem < J^ n b, (54) 
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where F is the fraction of the total particles emitted by PBH having value F < 0.2. Comparing above 
two equations, one can get 



0.38 



where r] evap — ^ h - is the baryon to photon ratio at evaporation. Assuming that rj is fixed from evapo- 
ration up to present times, i.e., r\ evav = ?7o and using the relation 770 ~ 2.8 x 10 _8 f2&/i 2 with flbh 2 fs 0.02, 
Eg. ([55]) becomes 



a evap < 1.064 x 10" 



The bounds imposed at nucleosynthesis can be converted into bounds on the initial PBH mass fraction 
for PBHs that evaporate at t evap — 400s: 

/ M- \ 1/2 

ai < 8 X 10 "° X (Wl) • (57) 
One may compare the constraint with that of standard cosmology (43l . |44| 

^<3xl0--(^) 1/2 . (58) 

Eauation (l5"Tj) further leads to 

3 „\V2 



on < 3.82 x 10" 19 



It turns out from the analysis by Carr et al. [25j that by considering the effects of quark and gluon 



emission by the PBHs, an earlier constraint [43J on the PBH mass fraction is actually weakened in the 
relevant mass range that we are also considering for our analysis in the Brans-Dicke case. 

B. Deuterium photodisintegration constraint 

The high-energy particles emitted by evaporating PBHs both during and after nucleosynthesis can 
be sufficiently energetic to disrupt primordial nuclei. One important reaction of this type is photo- 
disintegration which entails the destruction of primordial nuclei by high-energy PBH photons. Of all the 
primordial nuclei, deuterium is the most susceptible to photon-disintegration. If AM is the PBH mass 
evaporated between times t\ and £2 during which deuterons are destroyed and Mb is the baryonic mass, 
then 0] 

AM 

where / 7 is the fraction of mass that decays into photons, e is the depletion factor, and (3 are constants. 
t\ and t2 are taken to be the end of nucleosynthesis and the onset of recombination respectively. We 
estimate AM by the PBH mass evaporated shortly after nucleosynthesis. It is usually justified to take 




AM 

mT 



Ppbh 



Pb 



(61) 



with t evap some time after nucleosynthesis, either when a narrow mass range of PBHs evaporates, or 
straight after nucleosynthesis for an extended mass spectrum. Equation (|60[) then leads to 
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Since pi, oc a 3 and p ra d oc a 4 , it follows that 



Pb 



Prad 



a(t e ) 



Pb 



Prad 



te a(te) 



(63) 



as p ro< j = w Pc/2 at the time of matter radiation equality t e — 10 11 s. The baryon density 

parameter at equality is related to the present one, as the matter density parameter at equality is given 
by n m (t e ) « 1/2. 



Using equation (|64|) in equation (f62j). one gets 



2 n m (fo) 



a, 



/ 2 ^fcfa) 



(64) 



(65) 



Using numerical values of constants [26( / 7 = 0.1, E* = 10 1 Gev, (3 = 1, e ~ 1 and taking the ratio of 
the present baryonic density to the total matter density as 0.1, we get 



< 1.066 x 10 - 28 f^V /2 . 



For te-uap = 400 sec, the constraint becomes 



evap 



< 6.74 x 10" 



(66) 



(67) 



which gives 

a, < 5.1 x 10 
In standard case this constraint is [43, [451 



3/2 



«i < 10" 



(69) 



The recent analysis by Carr et al. [25j reveals that the standard deuterium constraint [43| is tightened 
by two orders of magnitude, since hadrodissociation of helium due to injected nuclei produces more 
deuterium. The same is expected to be true for the Brans-Dicke scenario, but, as in the earlier cases, 
only detailed numerical simulations could reveal the actual quantitative changes on the constraints. 



VII. SUMMARY AND CONCLUSIONS 



In this paper we have analyzed several astrophysical constraints on the initial mass fraction of primordial 
black holes that could grow due to accretion of radiation in the early Universe in Brans-Dicke cosmology. 
It is well-known that PBHs loose mass due to Hawking evaporation with their lifetime depending upon 
their initial masses. However, the accretion of radiation could be an efficient process in the BD formalism, 
depending upon the accretion efficiency, leading to the growth of mass and thereby enhancing their lifetime 
[THHl]]. Such PBHs, once formed, will influence the later cosmological epochs through their total density, 
and also through the products of their evaporation. Since the standard cosmological scenario is based 
on a sound observational footing, at least from the era of nucleosynthesis onwards, any modification 
due to surviving PBHs at various eras is tightly constrained by various observational results, such as 
the photon spectrum, CMB radiation, light element abundances, etc. Such constraints, in turn, can 
be translated into constraints on the initial PBH mass fraction. This translation is contingent on the 
particular cosmological evolution, and is thus sensitive to the theory of gravity considered, as has been 
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TABLE I: The variation of the upper bound of the initial PBH mass fraction(ai) with the accretion efficiency / 
for different cases is shown in the table 



Cause of the Constraint 


/ = o 




/ = 0.25 


/ = 0.45 


/ = 0.65 


Present Density 


3.43 x 10" 


■18 


1.69 x 10" 18 


0.62 x 10" 18 


0.01 x 10" 18 


Photon Spectrum 


5.34 x 10" 


26 


3.34 x 10 -26 


1.73 x 10 -26 


1.33 x 10 -27 


Distortion of CMB 


1.28 x 10" 


-21 


0.63 x 10" 21 


0.23 x 10" 21 


0.05 x 10" 22 


Helium abundance 


3.82 x 10" 


19 


3.01 x 10 -19 


2.18 x 10" 19 


0.60 x 10 -19 


Deuterium abundance 


5.10 x 10" 


21 


2.52 x 10~ 21 


0.94 x 10" 21 


0.02 x 10" 21 



shown earlier in the context of BD gravity without accretion [2(| , and also in the context of braneworld 
gravity [26| . Moreover, inclusion of the effect of accretion also impacts upon the constraints, as seen 
earlier in the context of braneworld gravity [26| , and also as shown by us in the present work. 

The summary of our results are presented in Table 1. Here we enlist the upper bounds on the initial 
mass fraction of PBHs that are allowed by taking into considerations various observational features as 
listed. These constraints obviously depend upon the accretion efficiency, as seen from the displayed 
numbers. The initial constraints are usually the most severe for those black holes whose lifetimes are 
comparable with the cosmic time of the epoch at which the observational constraint is imposed, and our 
displayed results correspond to such cases. We consider three different values of the accretion efficiency, 
with / « 2/3 corresponding to the maximum allowed in BD theory. The / = case corresponds to the 
BD formalism without accretion. It is seen that in all cases the inclusion of accretion strengthens the 
constraint on the upper bound of the initial PBH mass fraction. 

Since the present observable Universe is flat, the mass density of presently surviving PBHs should 
not exceed that of dark matter which is about 0.3 times the critical density. This imposes a constraint 
on initial PBH mass fraction which is of order 10 -18 for lower accretion efficiencies and grows by two 
orders for higher accretion efficiencies. The standard cosmology constraint is of order 10~ 18 . PBHs 
which evaporate between photon decoupling and the present age leave behind a spectrum that peaks at 
a temperature of the order of the black hole temperature at the onset of evaporation of PBHs with 
tevap ~ to- In this case the bound on the initial PBH mass fraction is comparable with standard 
cosmology for higher accretion efficiencies and grows weaker by one order for lower accretion efficiencies. 
If evaporation products are released around the Sunyayev-Zel'dovich time t$z ~ 4 x 10 -10 £o, they will 
fail to fully thermalize the background radiation. This time, however, is sufficiently early for the excess 
energy to distort the background blackbody spectrum. Limits on the allowed distortion of the CMB 
spectrum then imply limits on PBH mass fraction. For lower accretion efficiencies this constraint is of 
the same order as in standard cosmology, but increases more than three orders for maximally efficient 
accretion. If there were a population of PBHs evaporating during or after the era of nucleosynthesis, 
this could have led to significant change in the final light element abundances. Considering the Helium 
abundance as an example, we have found that the constraint on the initial PBH mass fraction once again 
varies as it started from a value of nearly one order less compared to standard cosmology and grows 
closer to it as accretion efficiency increases. Considering photon-disintegration and change in deuterium 
abundance, we find that the initial constraint is nearly of the same order as the standard value for low 
accretion efficiencies, increasing by nearly two orders for higher values. Comparing the constraints due to 
the different observational features considered in this work, we note that the photon spectrum imposes 
the most stringent limits on the initial PBH mass fraction. 

We conclude by emphasizing that Brans-Dicke cosmology which provides a viable alternative to the 
standard scenario, imposes upper bounds on the allowed initial mass fraction of primordial black holes, 
that are modified compared to standard cosmology. Depending upon the particular observed physical 
process used to impose the constraints, these upper bounds in BD gravity could either be strengthened or 
weakened compared to the case of standard gravity [2(| ■ It needs to be mentioned here that some of the 
constraints of standard gravity could themselves be modified by considering effects of quark and gluon 
emission and the resultant emission of secondary photons by PBHs, as discussed in the recent work by 
Carr et al. (25[. However, as shown in the present paper the inclusion of the effect of accretion tightens 
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the constraints in all cases since PBHs in BD gravity could start with a lower value of initial mass, and 
subsequently grow in size sufficiently [2l| to impact the observational features in future eras. Finally, it 
remains to be seen how additional effects such as considering PBHs in further general scalar-tensor models 
of gravity such as in [l8| , or taking into account the effects of backreaction of the PBHs on cosmological 
evolution, could modify the observational constraints on the initial mass spectrum. 
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